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ABSTRACT 

We present high resolution optical echelle spectra and IUE observations during a strong 
flare on 1993 December 22 in the very active, young, rapidly rotating, single K2 dwarf 
LQ Hya. The initial impulsive phase of the flare, which started sometime between 
2:42 UT and 4:07 UT, was characterized by strong optical continuum enhancement 
and blue-shifted emission lines with broad wings. The optical chromospheric lines 
reached their maximum intensity at «5:31 UT, by which time the blue-shift vanished 
and the optical continuum enhancement had sharply decreased. Thereafter, the line 
emission slowly decreased and the lines red-shift in a gradual phase that lasted at least 
two more hours. The Mg n lines behaved similarly. Quiescent C iv flux levels were 
not recovered until 21 hours later, though a data gap and a possible second flare make 
the interpretation uncertain. In addition to the typically flare-enhanced emission lines 
(e.g., Ha and H/3), we observe He I D 3 going into emission, plus excess emission (after 
subtraction of the quiescent spectrum) in other He I and several strong neutral metal 
lines (e.g., Mg I b). Flare enhancement of the far UV continuum generally agrees with 
a Si I recombination model. We estimate the total flare energy, and discuss the broad 
components, asymmetries, and Doppler shifts seen in some of the emission lines. 

Key words: stars: flare - stars: individual: LQ Hya - stars: activity - stars: chro- 
mospheres - stars: late-type - line: profiles 



1 INTRODUCTION 

LQ Hya (HD 82558) is a rapidly rotating (v sini = 25 km 
s _1 ), single K2 dwarf, classified as a BY Dra variable (Fekel 
et al. 1986a; Fekel, Moffett & Henry 1986b; Strassmeier & 
Hall 1988), with a photometric rotational period of 1.600881 
days (Strassmeier et al. 1997). Its high lithium abundance 
(Fekel et al. 1986a) suggests LQ Hya has age t < 7.5 x 10 7 
years (at least as young as the youngest Pleiades star); 
Vilhu, Gustafsson & Walter (1991) even suggest that it may 
be a pre-main sequence object. Saar, Piskunov & Tuominen 
(1992, 1994), Strassmeier et al. (1993) and Rice & Strass- 
meier (1998) found variable spot distribution in this star us- 
ing Doppler-imaging, and widespread magnetic fields have 



* Based on observations made with the William Herschel oper- 
ated on the island of La Palma by the Royal Greenwich Obser- 
vatory at the Spanish Observatorio del Roque de Los Muchachos 
of the Instituto de Astrofisica de Canarias, and on data obtained 
with the International Ultraviolet Explorer 



been detected by Saar et al. (1992, 1994), Basri & Marcy 
(1994), Donati et al. (1997), and Donati (1998). 

LQ Hya is also (not surprisingly) a very active star, 
as indicated by emission in several chromospheric lines, in- 
cluding Ca II H & K (Fekel et al. 1986a,b; Strassmeier 
et al. 1990), Ca II A8542 (Basri & Marcy 1994), and Ha 
(Vilhu et al. 1991), which can also appear as partly (Strass- 
meier et al. 1993) or completely filled-in absorption (Fekel 
et al. 1986a) . A filled-in He I D3 line is reported by Vilhu et 
al. (1991) and Saar et al. (1997). Strong UV chromospheric 
and transition region emission lines have also been found by 
Simon & Fekel (1987), and the star has been detected by 
ROSAT (Pye et al. 1995) and EUVE (Bowyer et al. 1996). 

Flares are believed to result from the release of magnetic 
free energy stored in the corona through reconnection (see 
reviews by Mirzoyan 1984; Haisch, Strong & Rodono 1991). 
Many types of cool stars flare (Pettersen 1989), sometimes 
at levels several orders of magnitude more energetic than 
their solar counterparts. In the dMe stars (or UV Cet type 
stars) optical flares are a common phenomenon, however, in 
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Table 1. Observing log WHT/UES (1993 December 22) 



UT JD ip S/N Description 

(h:m:s) (2449343.0+) 



01:14:26 


0.55 


0.444 


52 


Quiescent 


02:41:52 


0.61 


0.482 


108 


Quiescent 


04:35:48 


0.69 


0.532 


145 


Impulsive 


05:31:18 


0.73 


0.557 


118 


Flare maximum 


06:00:43 


0.75 


0.569 


125 


Gradual 


06:07:52 


0.76 


0.573 


134 




06:14:20 


0.76 


0.576 


132 




06:21:04 


0.76 


0.578 


129 




06:29:06 


0.77 


0.582 


123 




06:35:06 


0.77 


0.585 


124 




06:41:03 


0.78 


0.587 


129 




06:47:01 


0.78 


0.590 


130 




06:53:01 


0.79 


0.592 


122 




06:59:08 


0.79 


0.595 


123 




07:05:06 


0.80 


0.598 


122 




07:11:04 


0.80 


0.600 


128 




07:17:08 


0.80 


0.603 


128 




07:23:07 


0.81 


0.605 


124 




07:29:05 


0.81 


0.608 


129 





more luminous stars flares are usually only detected through 
UV or X-ray observations (e.g., Landini et al. 1986; Hiinsch 
& Reimers 1995; Ayres et al. 1994); optical flares are rare 
(Catalano 1990; Saar, Nordstrom & Andersen 1990; Henry 
& Newsom 1996). 

Ambruster & Fekel (1990) detected a strong ultraviolet 
flare on LQ Hya in four continuum bands between 1250 and 
1850 A, while no enhancement was evident in any of the 
chromospheric lines. Recently, HST GHRS observations by 
Saar & Bookbinder (1998) showed that many low-level flares 
are present in the transition region lines of this star. But 
consistent with its K2 spectral type and correspondingly 
increased optical continuum, strong optical flares are rare 
on LQ Hya (Montes et al. 1998b); Henry & Newsom (1996), 
for example, saw none. 

In this paper, we report one of these rare events: the 
detection of an unusually strong optical flare in LQ Hya 
through simultaneous observations of several optical chro- 
mospheric activity indicators: Ha, H/3, Na I Di, D2, He I 
D3, Mg 1 b triplet lines, and several UV chromospheric and 
transition region lines. 

In Sect. 2 we give the details of our observations and 
data reduction. In Sect. 3 we describe the different aspects of 
the flare deduced from our echelle spectra, such as the con- 
tinuum enhancement, the response of chromospheric lines 
to the flare, the variation of other photospheric lines, the 
energy released of various emission features as a function 
of time during the flare, and the line asymmetries. Finally, 
Sect. 4 gives the conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

Echelle spectroscopic observations of LQ Hya were obtained 
with the 4.2m William Herschel Telescope (WHT) and the 
Utrecht Echelle Spectrograph (UES) on 1993 December 22, 
covering several optical chromospheric activity indicators. 
These WHT/UES spectra were obtained with echelle 31 
(31.6 grooves per mm) and a 1024 x 1024 pixel TEK2 CCD 
as detector. The central wavelength is 5808 A covering a 



Table 2. Continuum variation* during the flare 



Obs. 


A4866 


A5175 


A5868 


A 6540 


(UT) 


(%) 


(%) 


(96) 


(%) 


02:42! 










04:36 


34 


32 


26 


23 


05:31 2 


10 


8 


7 


8 


06:01 


5 


5 


4 


6 


06:29 


6 


7 


5 


4 


07:29 


9 


9 


7 


2 



* near H/3, Mg I b, He I D3, and Ha, respectively 
Quiescent spectrum; 2 Flare maximum 



wavelength range from 4842 to 7725 A in a total of 44 
echelle orders. The reciprocal dispersion ranged from 0.048 
to 0.076 A/pixel. In Table @ we give the observing log. For 
each echelle spectrum we list the universal time (UT), the 
Julian Date (JD), the rotational phase (ip), and signal to 
noise ratio (S/N) obtained in the Ha line region. The ro- 
tational phase (ip) was calculated with the ephemeris re- 
cently given by Strassmeier et al. (1997) (To = 2445275.0, 
Pphtm — 1.600881), for alternative period determinations see 
also Strassmeier et al. (1993) and Jetsu (1993). 

The spectra have been extracted using the standard 
IRAF reduction procedures [J] (bias subtraction, flat-field di- 
vision, and optimal extraction of the spectra). The wave- 
length calibration was obtained using spectra of Th-Ar 
lamp. Finally, the spectra were normalized by a polynomial 
fit to the observed continuum. 

Frequent IUE SWP and LWP spectra of LQ Hya were 
taken between 15 and 24 Dec 1993. The data were reduced 
with standard IUE software and fluxes above background 
determined by simple integration (Table ^|). Here we ana- 
lyze the NEWSIPS calibrated data, resulting in some small 
changes in the measured fluxes relative to our initial results 
(Montes et al. 1998b). 



3 DESCRIPTION OF THE FLARE 

We detected a strong flare during the echelle observations 
of LQ Hya on 1993 December 22. The temporal evolution of 
the flare consists of an initial impulsive phase which started 
between 2:42 UT (last quiescent optical spectrum) and 4:07 
UT (end of the first IUE exposure with enhanced emission) . 
By the next optical spectrum (4:36 UT) strong increases in 
the chromospheric lines and continuum are seen. The opti- 
cal chromospheric lines reached maximum intensity at 5:31 
UT, by which time the continuum enhancement had already 
strongly decreased. After this, the lines slowly decreased in 
a gradual phase that lasted at least until the end of the 
observation (07:29 UT), i.e., > 2 hours. 

In the following we describe in detail the various as- 
pects of the flare deduced from our spectra, first exploring 
the time variation of the continuum, and then the response 

t IRAF is distributed by the National Optical Observatory, which 
is operated by the Association of Universities for Research in 
Astronomy, Inc., under contract with the NSF. 
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of the lines. Line emission is seen both in the "typical" chro- 
mospheric diagnostics (H I Balmer, and He I D3) and, after 
the subtraction of the quiescent spectrum, also in other He I 
lines (A4921.9, 5015.7, and 6678.2 A) and other strong lines 
such as the Na I Di and D2, the Mg I b triplet and several 
Fe 1 and Fe II lines. Finally, we calculate the energy release 
during the flare and we analyse the broad component and 
the asymmetry exhibited by some of the emission lines. 



3.1 The variation of the continuum 

Our echelle spectra show a change in depth of all the pho- 
tospheric lines due to continuum enhancement during the 
flare. We have determined the contribution of the flare to the 
total observed continuum by calculating what fraction of a 
linear continuum must be added to the quiescent spectrum 
in order to reproduce the corresponding flare spectrum. In 
Table ^ we give this contribution in representative spectral 
orders for several spectra during the flare. The maximum 
contribution of the flare to the continuum is reached at the 
beginning of the event (the impulsive phase) and decreases 
thereafter. It clearly depends on wavelength: in the H/3 line 
region the maximum contribution is 34 %, in the He I D3 
line region is 26 %, and in the Ha line region is 23 %, giv- 
ing a power law index F CO nt oc A~ 1,35 and an approximate 
blackbody temperature of T ~ 7500 K. 

The continuum behaviour is thus in agreement with 
photometry of other stellar flares, showing that the flare is 
initially dominated by strong continuum radiation, strongest 
at short wavelengths. The temperature indicated suggests 
that the plasma has already cooled somewhat and thus our 
"impulsive" spectrum (4:36 UT) may come somewhat late 
in the impulsive phase. 



3.2 The response of the optical Chromospheric 
lines to the flare 

To analyse the behavior of the various optical chromospheric 
lines during the flare, we first subtracted the quiescent spec- 
trum (UT:02:42), as is normally done in the analysis of flare 
stars. However, since the star is active, this procedure un- 
derestimates the total chromospheric contribution in these 
features, and ignores any variation (due to e.g., rotational 
modulation or evolution) of the "quiescent" state. To ob- 
tain total chromospheric contribution, we applied the spec- 
tral subtraction technique (i.e., subtraction of the rotation- 
ally broadened, radial- velocity shifted spectrum of a inactive 
star chosen to match the spectral type and luminosity class 
of LQ Hya; see Montes et al. 1995a, b, c). We used HD 10476 
(K1V) as the inactive reference star, taken from the spectral 
library of Montes & Martm (1998). 

In the case of the Ha and H/3 lines we have computed, 
for all the flare spectra, both the observed - quiescent (O-Q) 
and observed - reference (O-R) profiles (see Figures hi and 
^|). For the rest of the lines, not affected by chromospheric 
activity in the quiescent spectrum, we studied only the (O- 
Q) spectra (see Figures^, ^, |E| ^). In all cases, before the 
standard spectrum (quiescent or reference) is subtracted we 
take into account the contribution of the flare to the contin- 
uum (§3.1) in each spectral region, for each flare spectrum. 



3.2.1 The Ha line 

In Fig. ^ we have plotted the quiescent spectrum, the ref- 
erence star spectrum, and representative spectra during the 
development of the flare in the the Ha region. In the left 
panel we plot the observed spectra, in the central panel the 
(O-Q) profiles, and in the right panel the (O-R) profiles. This 
figure clearly shows the conspicuous Ha emission enhance- 
ment, from a weak emission above the continuum (quiescent 
spectrum) , to a very strong and broad emission at the max- 
imum of the flare. The excess Ha emission equivalent width 
in the (O-R) spectra, increases by a factor of ~2.7 in an 
interval of 2.8 hours. After the maximum the emission de- 
creases slowly; if modeled with an exponential decay EW oc 
exp(— 1//3), the e-folding time /3 ~2.5 hours in the first hour, 
slowing even further to /3 ~11 hours in the second hour. 

One significant aspect of these spectra is a broad emis- 
sion component in the both the (O-Q) and (O-R) profiles 
which makes the total emission poorly matched by a single- 
Gaussian fit. To study these profiles, we have therefore fitted 
them using two Gaussian components: a narrow component 
(N) having a FWHM of 56-69 km s _1 and a broad com- 
ponent (B) with 190 < FWHM < 293 km s -1 . In Table | 
we list the parameters (I, FWHM, EW) of the broad and 
narrow components. As can be seen in this table, the con- 
tribution of the B component to the total EW and FWHM 
of the line is a maximum in the impulsive phase. The line 
profiles are also asymmetric, and the two Gaussian fit is op- 
timized when the broad component is blue- or red-shifted 
with respect to the narrow component (see the AA = An - 
Ab value in Table ^|). We discuss this further in §3.5. 

In Table ^| we also give, for the total subtracted spec- 
tra, the peak emission intensity (I), the excess Ha emission 
equivalent width (EW(Ha)), and absolute fluxes at the stel- 
lar surface, logFs(Ha), in erg cm -2 s _1 obtained with the 
calibration of Hall (1996). The time evolution of the EW( 
Ha) during the flare is shown in Fig. ^| 

3.2.2 The H/3 line 

Fig. |^ is as Fig. |l| for the H/3 line. In this case the line 
changes from mostly filled-in absorption (quiescent), to a 
strong and broad emission line at flare maximum. The excess 
H/3 emission EW in (O-R) spectra increases by a factor of 
5.8 from the quiescent level to the maximum, a considerably 
larger enhancement than in the Ha line. However, during the 
gradual phase the emission declines slightly more rapidly 
than the Ha emission (see Fig. ^|), with f3 = 2.3 hours in the 
first hour of decay, slowing to /3 = 5.9 hours in the second. 
The more rapid decay of H/3 has also been observed in other 
solar and stellar flares (Johns-Krull et al. 1997). 

Like Ha, very broad wings are visible in the (O-R) H/3 
profiles, and the results of double Gaussian fits are given in 
Table ^j. The contribution and FWHM of the broad compo- 
nent again reach a maximum in the impulsive phase and the 
profiles show an asymmetry similar to that seen in Ha. 

3.2.3 The Ha/ H/3 Balmer decrement 

The Balmer decrement can be an important indication of 
physical parameters such as electron density (Kunkel 1970). 
It is commonly tabulated for stellar flares with respect to the 
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H7 line (Hawley & Pettersen 1991). Lacking H7 data in our 
spectra, we have instead calculated the Balmer decrement 
from the EW(Ha)/EW(H/3) ratio, assuming that the LQ 
Hya continua at Ha and H/3 have a ratio appropriate to a 
blackbody at T e// = 4900 K (F A6 « 3 /Fa486i = 1-08). The 
values obtained are given in Table m 

We found that the Ha/H/3 Balmer decrement is shal- 
lower during the flare, changing from 3.15 at the quiescent 
state to 1.46 at the impulsive and maximum phases of the 
flare, indicating (unsurprisingly) a significant change in the 
properties of hydrogen emitting regions during the flare. In 
all the flares on dMe stars tabulated by Hawley & Pettersen 
(1991), the decrement for H/3 and higher Balmer lines was 
also shallower during the flare. The Ha/H/3 ratio can fall 
to < 1 during large stellar flares (Kunkel 1970) and in the 
largest solar flares (Zirin & Ferland 1980). However, in other 



stellar flares steeper decrements (2.0) are reported (Worden 
et al. 1984) and in the solar flares observed by Johns-Krull 
et al. (1997) this ratio only changed from 2.5 to 1.45. 

3.2.4 The He I D3 and other He I lines 

In Fig. ^ we show the quiescent spectrum and representa- 
tive flare spectra in the the He I D3 region; the left panel 
displays the observed spectra and in the right panel, the 
(O-Q) profiles are plotted. The He I D3 line, completely 
filled-in in quiescence, goes into emission during the flare, 
reaching a maximum intensity (0.346) at the same time that 
the Balmer lines. The subtracted profiles again show broad 
wings; an analysis like that used with Ha and H/3 yields the 
parameters in Table [|. 

This line is a well known diagnostic of flares in the Sun 
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Figure 2. H/3 observed spectra (left panel), after the subtraction of the quiescent spectrum (central panel) and after the spectral 
subtraction (right panel) 



and in other stars. In the Sun the D3 line appears as absorp- 
tion in plage and weak flares and as emission in strong flares 
(Zirin 1988). The He I D3 feature is typically in emission in 
dMe stars (e.g., Pettersen, Evans & Coleman 1984); in flare 
like events in UV Ceti flare stars (Kunkel 1970; Bopp & 
Moffett 1973); in strong flares in more luminous, active stars 
such as the RS CVn systems II Peg (Huenemoerder & Ram- 
sey 1987; Montes et al. 1997; Berdyugina, Ilyin & Tuominen 
1998) and UX Ari (Montes et al. 1996b; 1997), the weak- 
lined T Tauri star V410 Tau (Welty & Ramsey 1998), and 
in the active G5V k Cet (Robinson & Bopp 1987). Thus, 
the detection of prominent D3 emission indicates that we 
are observing a very strong flare in LQ Hya. 

Other He I lines have been also reported in stellar flares 
(Bopp & Moffett 1973; Hawley & Pettersen (1991); Abdul- 
Aziz et al. 1995; Abranin et al. 1998). In our spectra, after 



the subtraction of the quiescent spectra, we have also found 
an excess emission in other He I lines at 4921.93, 5015.68 
and 6678.15 A. In particular, He 1 A6678.15 appears super- 
imposed on the Fe I A6677.99 absorption line and the excess 
emission can be seen even in the observed spectra (see Fig. ^ 
left panel) . When the quiescent spectra is subtracted, excess 
He 1 emission is clearly seen (see Fig. [| right panel), with 
a maximum peak intensity of 0.17. The profiles have been 
again fitted with two Gaussians; the corresponding param- 
eters are given in Table [j. The temporal evolution of this 
He 1 line during the flare is similar to the D3 line. 

3.2.5 The Na I Di and D 2 lines 

The Na I Di and D2 resonance lines at 5895.9 A and 5890.0 A 
are collisionally-controlled in the atmospheres of late-type 
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Table 3. Ha parameters during the flare in the subtracted profiles for the two Gaussian component fits and for the total emission 





Ho broad component 




Ha narrow component 










Obs. 


I 


FWHM 


EW B 


b/t 


I 


FWHM 


EW N 


n/t 


AA 


1 T 


EW T 


log Ft 


(UT) 




(A) 


(A) 


(%) 




(A) 


(A) 


(%) 


(An - Ab) 




(A) 




Obs. - Quiescent 


02:42 (Quiescent) 




















0.000 


0.000 


0.00 


04:36 


0.137 


7.709 


1.098 


78.4 


0.117 


2.414 


0.304 


45.2 


+0.34 


0.254 


1.401 


6.70 


05:31 (Maximum) 


0.192 


6.766 


1.371 


65.4 


0.281 


2.427 


0.725 


34.6 


-0.06 


0.471 


2.096 


6.88 


06:01 


0.174 


6.021 


1.113 


66.1 


0.240 


2.227 


0.570 


33.9 


-0.36 


0.417 


1.683 


6.78 


06:29 


0.140 


6.414 


0.948 


66.7 


0.199 


2.238 


0.473 


33.3 


-0.42 


0.331 


1.421 


6.71 


07:29 


0.121 


6.567 


0.837 


64.3 


0.193 


2.267 


0.465 


35.7 


-0.66 


0.315 


1.302 


6.67 


Obs. - Reference 


02:42 (Quiescent) 


0.053 


4.168 


0.236 


20.3 


0.714 


1.221 


0.927 


79.7 


-0.12 


0.753 


1.163 


6.62 


04:36 


0.194 


6.414 


1.316 


59.4 


0.615 


1.371 


0.898 


40.6 


+0.21 


0.802 


2.214 


6.90 


05:31 (Maximum) 


0.307 


5.532 


1.809 


58.5 


0.835 


1.442 


1.281 


41.5 


+0.07 


1.132 


3.090 


7.04 


06:01 


0.270 


5.153 


1.481 


54.4 


0.834 


1.397 


1.240 


45.6 


-0.11 


1.085 


2.721 


6.99 


06:29 


0.226 


5.326 


1.280 


51.2 


0.835 


1.221 


1.221 


48.8 


-0.17 


1.043 


2.501 


6.95 


07:29 


0.201 


5.449 


1.167 


48.6 


0.822 


1.397 


1.233 


51.4 


-0.09 


1.006 


2.400 


6.93 



Table 4. H/3 parameters during the flare in the subtracted profiles for the two Gaussian component fits and for the total emission 





H/3 broad component 




H/3 


narrow 


component 










Obs. 


I 


FWHM 


EW B 


b/t 


I 


FWHM 


EW N 


n/t 


AA 


It 


EW T 


Ha/H/3 


(UT) 




(A) 


(A) 


(%) 




(A) 


(A) 


(%) 


(An - Ab) 




(A) 




Obs. - Quiescent 


02:42 (Quiescent) 




















0.000 


0.000 


0.00 


04:36 


0.208 


5.185 


1.142 


86.3 


0.117 


1.450 


0.181 


13.7 


+0.03 


0.326 


1.323 


1.143 


05:31 (Maximum) 


0.279 


4.440 


1.142 


69.7 


0.357 


1.507 


0.572 


30.3 


-0.08 


0.644 


1.888 


1.498 


06:01 


0.203 


4.429 


0.957 


65.5 


0.328 


1.439 


0.503 


34.5 


-0.15 


0.532 


1.450 


1.253 


06:29 


0.178 


4.599 


0.870 


70.7 


0.248 


1.361 


0.359 


29.3 


-0.20 


0.429 


1.230 


1.246 


07:29 


0.158 


4.569 


0.767 


73.9 


0.178 


1.428 


0.271 


26.1 


-0.30 


0.366 


1.038 


1.353 


Obs. - Reference 


02:42 (Quiescent) 




















0.437 


0.398 


3.152 


04:36 


0.232 


5.102 


1.253 


76.5 


0.359 


1.004 


0.384 


23.5 


-0.09 


0.610 


1.636 


1.460 


05:31 (Maximum) 


0.327 


4.308 


1.501 


65.5 


0.663 


1.118 


0.789 


34.5 


-0.13 


1.003 


2.290 


1.455 


06:01 


0.250 


4.343 


1.156 


60.4 


0.661 


1.076 


0.757 


39.6 


-0.23 


0.920 


1.913 


1.534 


06:29 


0.211 


4.627 


1.036 


61.4 


0.594 


1.028 


0.650 


38.6 


-0.30 


0.814 


1.686 


1.600 


07:29 


0.191 


4.713 


0.955 


63.5 


0.514 


1.000 


0.548 


36.5 


-0.28 


0.718 


1.503 


1.723 



stars and thus provide information about chromospheric ac- 
tivity (see Montes et al. 1996b, 1997, 1998a). Recent models 
of these lines for M dwarfs have been made by Andretta, 
Doyle, & Byrne (1997). Using model solar flares, Falchi, Fal- 
ciani & Smaldone (1990) found that the Na I D2 line shows 
modifications only in the core of its profile. 

Subtraction of the modified reference star from the qui- 
escent spectrum reveals that excess emission in the Na I Di 
and D2 lines is very small; we have therefore analysed only 
the (O-Q) spectra to study the flare in these lines. As can be 
seen in Fig. pi the lines show a clear filling-in of the core its 
profile, reaching its maximum intensity (0.13) at the same 



time as the other chromospheric lines. The peak intensity 
(I), FWHM, and EW for both lines are given in Table §. 

3.2.6 The Mg I b triplet lines 

The strong Mg I b triplet AA5167, 5172, 5183 is formed in the 
lower chromosphere and the region of temperature minimum 
and they are good diagnostics of activity (Basri, Wilcots & 
Stout 1989; Gunn & Doyle 1997; Gunn, Doyle & Houdebine 
1997). In some stellar flares this lines exhibit a central re- 
versal (Bopp & Moffett 1973; Mochnacki & Schommer 1979; 
Abdul-Aziz et al. 1995; Abranin et al. (1998). In the qui- 
escent spectrum very strong absorption lines are observed, 
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Figure 3. He I D3 observed spectra (left panel) and after the 
subtraction of the quiescent spectrum (right panel) 
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6074 6676 



Wavelength (A) 



Observed - Quiescent 




Broad component 
Narrow component 



Wavelength (A) 



Figure 4. He I A6678 observed spectra (left panel) and after the 
subtraction of the quiescent spectrum (right panel) 



without evidence of filling-in by activity. In the flare spec- 
tra, however, a small reversal is observed in line cores (Fig. ^ 
left panel). After the subtraction of the quiescent spectra, 
excess emission is clearly visible, both in the three Mg I b 
lines (Fig. |^ right panel), and (even more intensely) in the 
Fe 11 feature at 5169.0 A. The measured parameters for these 
lines are given in Table ^. 

3.2.7 The variation of other photo spheric lines 

The upper photosphere of the late-type stars is also affected 
by the chromospheric activity, as seen in the filling-in of 
strong photospheric lines with relatively low excitation po- 
tentials in active dwarfs (Basri et al. 1989) and pre-main 
sequence stars (Finkenzeller & Basri 1987). During mod- 
erately strong solar flares a large number of photospheric 
lines are filled-in when the quiescent spectrum is subtracted 
(Acampa et al. 1982; Mauas 1990; Johns-Krull et al. 1997). 
Filling-in of some metal lines are also reported in earlier ob- 
servations of flares of the UV Cet-type stars by Mochnacki 
& Schommer (1979; and references therein), Hawley & Pet- 
tersen (1991), and in recently observations by Abdul-Aziz 
et al. (1995) and Abranin et al. (1998). 

We have observed a slight filling in of many photo- 
spheric absorption lines in our spectra during the flare. 
These include all the lines reported to be activity-sensitive 



by Basri et al. (1989) and Finkenzeller & Basri (1987), 
many of the lines appearing in the high resolution spec- 
tra of a moderate solar flare (Johns-Krull et al. 1997), all 
the lines identified in the spectra of a more energetic solar 
flare (Acampa et al. 1982), and some of the lines reported 
in stellar flares (Mochnacki & Schommer 1979; Abdul- Aziz 
et al. 1995; Abranin et al. (1998). The lines with the largest 
filling in are those of multiplet 42 of Fe 11 (4943.9, 5018.4, 
5169.0 A). Other lines with significant filling are the lines 
of multiplet 48 (5316.8, 5362.9 A) and 49 (5197.6, 5234.6, 
5276.0, 5316.6 A) of Fe 11. 



3.3 The flare in the UV 

Ultraviolet line and continuum fluxes are given in (Table ^) . 
We coadd selected strong chromospheric (CHR) and transi- 
tion region (TR) line fluxes to improve statistics: /tr is the 
sum of N v (T form ~ 1.25 x 10 5 K), Si IV (T form ~ 8 x 10 4 
K) and C IV (Tf orm ~ 10 s K), while /chr coadds O I and 
C 1 (T form ~ 7 x 10 3 K), C 11 (T form ~ 2 x 10 4 K), and 
Si 11 (T form ~lx 10 4 K). Figure m shows the quiescent (an 
average of 10 exposures), impulsive (JD sta rt = 2449343.568; 
1:37 UT), and just post-maximum (= "UV peak"; JD sta rt = 
2449343.783; 6:48 UT) UV spectra. Figure | depicts Mg 11 
spectra between the SWP data in Figure plus an average 
quiescent spectrum. Figure nC] shows the time variation of 
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Table 5. He I D3 parameters during the flare in the (Observed - Quiescent) profiles for the two Gaussian component fits and for the 
total emission 



He I D3 broad component He I D3 narrow component 



Obs. 


I 


PWHM 


EW B 


b/t 


I 


FWHM 


EW N 


N /t 


AA 


It 


EW T 


(UT) 




(A) 


(A) 


(%) 




(A) 


(A) 


(%) 


(A N - A B ) 




(A) 


02:42 (Quiescent) 




















0.000 


0.000 


04:36 


0.033 


3.400 


0.117 


43.0 


0.160 


0.912 


0.155 


57.0 


+0.32 


0.198 


0.272 


05:31 (Maximum) 


0.051 


2.981 


0.161 


46.5 


0.226 


0.772 


0.185 


53.5 


+0.03 


0.272 


0.346 


06:01 


0.049 


2.450 


0.128 


53.3 


0.142 


0.745 


0.112 


46.7 


-0.13 


0.190 


0.241 


06:29 


0.030 


2.802 


0.088 


43.1 


0.133 


0.825 


0.117 


56.9 


-0.21 


0.161 


0.205 


07:29 


0.034 


2.554 


0.093 


42.0 


0.145 


0.826 


0.128 


58.0 


-0.11 


0.184 


0.220 



Table 6. He I A6678 parameters during the flare in the (O - Q) profiles for the two Gaussian component fits and for the total emission 



He I A6678 broad component He I A6678 narrow component 



Obs. 


I 


FWHM 


EW B 


b/t 


I 


FWHM 


EW N 


n/t 


AA 


1 T 


EW T 


(UT) 




(A) 


(A) 


(%) 




(A) 


(A) 


(%) 


(An - A B ) 




(A) 


02:42 (Quiescent) 




















0.000 


0.000 


04:36 


0.094 


0.911 


0.091 


95.4 


0.017 


0.249 


0.004 


04.6 


-0.08 


0.106 


0.095 


05:31 (Maximum) 


0.093 


0.857 


0.085 


72.4 


0.079 


0.384 


0.032 


27.6 


-0.09 


0.171 


0.117 


06:01 


0.080 


0.797 


0.067 


84.5 


0.041 


0.285 


0.012 


15.5 


-0.08 


0.120 


0.080 


06:29 


0.062 


0.833 


0.055 


78.3 


0.040 


0.360 


0.015 


21.7 


-0.13 


0.101 


0.070 


07:29 


0.047 


1.060 


0.053 


68.6 


0.060 


0.381 


0.024 


31.4 


-0.18 


0.111 


0.078 



Table 7. Na I D line parameters during the flare in the (O - Q) profile 







Na 1 D 2 






Na 1 Di 




Obs. 


I 


FWHM 


EW 


I 


FWHM 


EW 


(UT) 




(A) 


(A) 




(A) 


(A) 


02:42! 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


04:36 


0.058 


0.674 


0.036 


0.078 


0.549 


0.038 


05:31 2 


0.106 


0.584 


0.059 


0.127 


0.487 


0.065 


06:01 


0.095 


0.600 


0.057 


0.104 


0.528 


0.056 


06:29 


0.074 


0.718 


0.048 


0.092 


0.657 


0.052 


07:29 


0.090 


0.915 


0.063 


0.117 


0.666 


0.074 



Quiescent spectrum; 2 Flare maximum 



Table 8. Mg I b triplet and Fe II A5169 line parameters during the flare in the (O - Q) profile 







Mg 1 b 3 






Mg 1 b 2 






Mg 1 bi 






Fe 11 A5169 




Obs. 


I 


FWHM 


EW 


I 


FWHM 


EW 


I 


FWHM 


EW 


I 


FWHM 


EW 


(UT) 




(A) 


(A) 




(A) 


(A) 




(A) 


(A) 




(A) 


(A) 


02:42! 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


04:36 


0.078 


0.531 


0.044 


0.069 


0.503 


0.037 


0.060 


0.467 


0.030 


0.154 


0.458 


0.075 


05:31 2 


0.137 


0.392 


0.057 


0.122 


0.379 


0.049 


0.111 


0.350 


0.041 


0.275 


0.401 


0.118 


06:01 


0.111 


0.425 


0.050 


0.102 


0.381 


0.042 


0.089 


0.359 


0.034 


0.222 


0.402 


0.095 


06:29 


0.093 


0.342 


0.034 


0.082 


0.349 


0.030 


0.069 


0.337 


0.025 


0.177 


0.416 


0.078 


07:29 


0.103 


0.414 


0.045 


0.100 


0.318 


0.034 


0.082 


0.265 


0.023 


0.204 


0.387 


0.084 



Quiescent spectrum; 2 Flare maximum 
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Figure 5. Spectra near the Na I Di and D2 lines before (left 
panel) and after subtracting the quiescent spectrum (right panel) 




5165 5170 5175 51B0 51S5 5165 5170 5175 5150 5185 



Wavelength (A) Wavelength (A) 

Figure 6. Spectra near the Mg I b triplet lines before (left panel) 
and after the subtraction of the quiescent spectrum (right panel) 



/tr and /chr during the period of the optical flare. We note 
that while the UV "impulsive" spectrum ends before the first 
flare-affected optical spectrum (4:36 UT), it already shows 
noticeably enhanced (1. 5 x quiescent) transition region (TR) 
and continuum emission. Thus the true "impulsive" phase 
starts earlier than the optical data indicates, with high Tf orm 
TR lines leading the evolution (as is typical). 

Besides the lines measured in Table [], numerous weak 
lines are also present (typically in UV peak and/or quies- 
cent spectra). Most certain of these are C I 1261, 1277, 1280, 
1560 A, Si in A 1295-98, 1892 A, S I 1820 A, Al in 1855, 
1863 A, Fe ii 1611, 1633, 1671 A, and O in 1666 A. The 
following are also detected with less certainty: C I 1463, 
1493 A, S i 1473, 1483, 1487, 1900 A, Si I 1258, 1267, 1565, 
1574 A, Si ii 1533A, Fe n 1372, 1588, 1598 A, S iv 1417 A, 
O IV 1407 A, O v 1371 A, and S v 1502 A. Also of interest 
is the possible detection of the coronal 1354.1 A Fe xxi line 
(Tform ~ 10 7 K) in the UV peak and (perhaps) in the qui- 
escent spectrum. Possible contributions from other nearby 
lines (e.g., C I 1354.3 A) make this identification uncertain, 
though. 

The strongest high Tf orm line, C iv, shows a strong blue- 
ward asymmetry in the impulsive spectrum (its centroid 
is shifted by ~ —250 km s" 1 ), perhaps foreshadowing the 
blueshifts in the broad components of the optical Balmer 
and He I lines seen later (see §3.5). Si IV and N V (with sim- 



ilar Tf orm ) also appear to be blue-shifted, though the shift 
is uncertain due to probable blends and the weakness these 
features. Cooler strong lines (C II, He II) show no shift. The 
lack of H2 rules out the 1547 A line as responsible, but we 
note that numerous features near C IV are coincident with 
(normally weak) Si I (1545, 1552, 1555-8, 1563 A), and C I 
(1542 A) lines. Thus, the apparent C IV blueshift may be 
partly due to blends. 

There are no SWP data during optical flare maximum, 
though two LWP spectra appear to show Mg 11 peaking at 
about the this time (Fig. H). The presence of ISM absorp- 
tion, emission self-reversal, and nearby blends complicates 
the interpretation of these features. We first constructed a 
"quiescent" spectrum Fq from the average of eight LWP 
images which appeared uncontaminated by flares. We then 
removed the ISM absorption from Fq by adding a 2 pixel 
Gaussian at the center of the absorption feature. We mod- 
eled the Mg 11 lines in the flare very simply, with Ffl aro = 
AFq + BG B - CGism + D, where Fq is the ISM-corrected 
quiescent spectrum, Gb and Gism are broad and narrow (2 
pixel) Gaussians, and A, B, C , and D are adjustable con- 
stants. The width and central A of Fq and Gb were also 
allowed to vary; results are shown in Fig. ^. The advantage 
of using Fq as a template is that we account for the intrinsic 
Mg II shape and (partially) remove blends. 

We find that the first Mg 11 spectrum shows a 2.5 x 
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Figure 7. IUE SWP-lo spectra showing fluxes at earth (per A) for the mean quiescent phase (average of 10 spectra; heavy solid), 
the early impulsive phase (JD 2449343.568; thin solid) and the near-maximum phase (JD 2449343.783; dotted). The spectra have been 
smoothed by a 3 pixel running mean. 




Wavelength (A) 



Figure 8. IUE LWP-hi spectra of the Mg II k line showing fluxes 
at earth (per A) for the mean quiescent phase (average of 8 spec- 
tra; solid) , and two spectra near flare peak (taken between the two 
SWP spectra in Fig. Q); the later spectrum is offset by 10" 11 for 
clarity. Spectra with the ISM absorption approximately removed 
(dashed) and simple models combining ISM, enhanced quiescent, 
and broad components (see text) are shown (dotted); blends are 
also noted (thick vertical marks), as is the broad component peaks 
(note the shift between the spectra). 



enhancement over quiescent in its narrow component, and 
a broad («250 km s _1 FWHM) component with a flux of 
/b = 3.3 x 10" 12 ergs cm -2 s" 1 at earth in Mg II k. The 
broad component in this spectrum is Doppler shifted by 
~ —40 km s _1 . The continuum is also significantly enhanced 



(by at least lOx), with (/con) ~ 2 x 10~ 12 ergs cm" 2 s" 1 
A" 1 at 2800 A. The second spectrum, taken about 1.5 hours 
later, showed both broad and narrow components reduced 
by a factor of w 1.5, with the broad component now ba- 
sically unshifted but with a similar width (Fig. []). Thus, 
the Mg II lines respond to the flare much like the optical 
chromospheric features, initially exhibiting enhanced emis- 
sion with broad, blue-shifted components, which gradually 
drift to the red and weaken as the flare evolves. 

The next SWP spectrum (UV peak, coinciding the with 
early gradual phase in the optical data) reveals > 20 x en- 
hancements in the TR and > 4x enhancements in the CHR. 
Several high Tf orm lines only weakly detected in the quies- 
cent spectrum (C III 1175 A, N v 1240 A) are also greatly 
strengthened. Several lines not listed in Table ^| appear in 
this spectrum; the more certain of these include the C I com- 
plexes (1261 A, 1277 A, and 1280 A) and the Si ill complex 
at 1295-1299 A, (combined fluxes / « 1.3 and 0.9 x 10" 13 
ergs cm -2 s _1 ), O in 1666 A, Fe n 1671 A, Al in 1854+1862 
A, and Si ill 1892 A (with / ta 0.8, 1.3, 2.4, and 1.3 x 10" 13 
ergs cm -2 s _1 , respectively). 

The UV continuum is also noticably enhanced in this 
spectrum. Machado & Henoux (1982) suggested similar en- 
hancements in solar flares at A < 1524 A and A < 1682 Aare 
due to are due bound-free radiation from the ground and 
first excited levels of Si I excited by UV line emission (pri- 
marily by C II and C IV, respectively). Exploring this idea, 
Phillips, Bromage & Doyle (1992) studied the power per 
A in 50 A (relatively line free) bands centered at 1490 A and 
1595 A in stellar flares, and found that they correlated well 
with line power in C II and C IV, respectively: specifically, 
logP C ont(1490A) = 0.94 logPc ii + 0.3 and log P cont (1595A) 
= 1.04 logPc iv — 2.9. We have integrated the flux in sim- 
ilar bands (/ CO nt in Table note that the small / con t in 
the non-flare spectra may be due to the sum of many weak 
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Figure 9. The change of EW of several optical chromosphcric 
lines during the flare 



lines. We find good agreement for the 1490 A continuum, 
with log P CO nt( 1490 A) (= /cont/50A xAnd 2 , where d is the 
stellar distance) within an average 0.10 dex of the predic- 
tion for the impulsive and UV max spectra. Agreement is 
less good for log P C ont (1595 A), which is overestimated by 
an average of 0.31 dex. We see no Si I edges at 1524 A or 
1682 A, but this is perhaps not surprising given the low res- 
olution of IUE (Phillips et al. 1992). In general, though, the 
far UV flare continuum of LQ Hya seems consistent with the 
Si I recombination model. 

If LQ Hya were very young and still had significant 
circumstellar material, one might expect X-ray excitation 
of the gas in such a strong flare. On the suggestion of the 
referee, we explored whether any of the weak features in the 
impulsive or UV peak spectrum might be due to excitation 
of circumstellar H2. As there is no evidence in these spectra 
for the four strongest H2 features seen in T Tau (1446, 1490, 
1505, 1562 A; Brown, de M. Ferraz & Jordan 1984), it seems 
unlikely there is much nearby gas. This makes it less likely 
that LQ Hya is pre-main sequence, as suggested by Vilhu et 
al. (1991). 

The UV fluxes do not return to their quiescent level 
until almost a day later (JD 2449344.744), though lack of 
data for ~0.8 day after the UV maximum makes the inter- 
pretation ambiguous. The time between UV maximum and 
quiescence covers over half a rotation, making it difficult for 
the flare (if localised) to have remained visible unless it was 
very near the pole, or in a very extended (> i?*) loop. The 
sharp drop in flux from 344.589 to 344.744 could then be 
due in part to the flare's finally disappearing over the limb. 
Perhaps more likely, the enhancement seen on JD 2449344 
might be due to a second flare and the fast decay a conse- 
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Figure 10. Evolution of the combined IUE chromosphcric fluxes 
(at earth) of O 1 (1304 A), C n (1335A), C 1 (1657 A), and Si 11 
(1810 A) (=/ CHR , o), chromospheric Mg n/5 (2800 A; *), the 
combined transition region fluxes of N v (1240 A), Si iv (1400A) 
and C IV (1550A) (=/tRi a ) 5 and the UV continuum flux (aver- 
age of the total / CO nt in two 50 A bands centered at 1490 A and 
1595 A; A). Horizontal solid lines through the data points indi- 
cate the exposure durations. The first time points (plotted arbi- 
trarily at 0.4) show the mean quiescent fluxes and their errors 
(vertical line). The span of the optical observations is also indi- 
cated, with the line type indicating the spectrum's appearance 
(dashed=quiescent, heavy solid=impulsive/flare peak, thin solid 
= gradual); the optical maximum is marked with a tick. 



3.4 Estimate of energy released 

To estimate the flare energy released in the observed optical 
chromospheric lines, we converted the EW into absolute sur- 
face fluxes and luminosities. Since we have not observed the 
entire flare, and are missing important lines (e.g., the satu- 
rated Ly a; Mg II; He I 10830 A) our estimates are only lower 
limits to the total flare energy in chromospheric lines. We 
have used the calibration of Hall (1996) to obtain the stellar 
continuum flux in the Ha region as a function of (B - V) and 
then convert the EW into absolute surface flux. For the other 
lines, we have used the continuum flux at Ha corrected for 
the difference in the continuum Fa6563/Fa, assuming Fa is 
given by a blackbody at T e ff = 4900 K. The corresponding 
absolute fluxes at flare maximum (erg cm -2 s _1 ), and total 
flux (ergs cm -2 ) integrated over the the observation interval 
(~ 3 h) are given in Table [lo| We converted these fluxes into 
luminosities using the radius R = 0.76 R0 (Strassmeier et 
al. 1993). We have not estimated the energy released in the 
UV chromospheric lines, since, without the (saturated) Lya 
and Mg 11 lines, the total /chr will be greatly underesti- 
mated. We find that the estimated total energy released in 



quence of its lower energy (E a3ble oc t decay ; Lee, Petrosian & optical chromospheric lines is J5chr > 5.7 10 33 ergs, indicat- 



McTiernan 1993; Saar & Bookbinder 1998). 



ing that this LQ Hya flare is more energetic in line emission 
than an average flare on a dMe star, where typically 10 28 
ergs < -Echr < 10 34 ergs (see Hawley & Pettersen 1991). 
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Table 9. UV continuum fluxes, chromospheric (CHR) and transition region (TR) lines fluxes 



JD s tart 


tcxp 








/ (io- 


13 ergs 


cm 2 s 


at earth) 1 






/CHR 


/tr 


,/cont 


/cont 


-2449000 


00 


C in 


N v 


O i 


C II 


Si iv 


C iv 


tic II 


f~* T 
V I 


Si II 


Mg II 2 


sum 


sum 


(1490A) 


(1595A) 


342.742 


2400 
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3.5 


2.0 
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1.0: 
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1 lines at 1175A, 1239+43A, 1302+5+6A, 1335A, 1394+1403A, 1548+51A, 164lA, 1657A, 1808+17A, 2796+2803A, respectively. 

2 uncorrected/corrected for ISM absorption 3 not shown in Fig. 10. 4 colon = uncertain measurement 5 noisy spectrum 
6 line saturated, / > 38.0 X 10~ 13 ergs cm -2 s — 1 , flux estimated from Gaussian fit to line wings and shoulder. 



Table 10. The flare energy released the chromospheric lines 



Line 


F^(max) 
(106) 
(ergs cm -2 s — 1 ) 


jF x dt 
(10 10 ) 

(ergs cm -2 ) 


(max) 

(10 29 ) 

(org s -1 ) 


jh x dt 
(10 33 ) 

(erg) 


H/3 


6.269 


5.876 


2.204 


2.066 


Mg I b 3 


0.200 


0.210 


0.070 


0.074 


Fe ii A5169 


0.414 


0.406 


0.145 


0.143 


Mg I b 2 


0.172 


0.174 


0.062 


0.061 


Mg I bi 


0.144 


0.137 


0.051 


0.048 


He I D 3 


1.270 


1.234 


0.447 


0.434 


Na I D 2 


0.217 


0.253 


0.076 


0.089 


Na I Di 


0.239 


0.278 


0.084 


0.098 


Ha 


7.510 


7.242 


2.641 


2.546 


He I A6678 


0.415 


0.412 


0.146 


0.145 


Total lines 


16.85 


16.23 


5.926 


5.704 



3.5 Line asymmetry 

Very broad wings have been found in Mg II and in the sub- 
tracted profiles of Ha, H/3, He I D3 and He I A6678 lines, and 
double Gaussian (narrow=N and broad=B) fits were used 
to model them. The contribution of the broad component 
to the total EW and FWHM of these lines reaches a max- 
imum in the impulsive phase and then decreases. The line 
profiles are also asymmetric, with the B component often 
shifted relative to the N. In the impulsive phase we found 
a blue asymmetry in the lines - the B component appears 
blue-shifted with respect to the N. At the maximum of the 
flare AA « 0, and during the gradual phase a red asymmetry 
is present, with AA increasing with time (see Tables 3-6). 

Similar broad Balmer emission wings have been seen in 
dMe stars and chromospherically active binaries connected 
with flare events, with similar blue or red asymmetries in 
some cases (Eason et al. 1992; Phillips et al. 1988; Doyle et 
al. 1988; Gunn et al. 1994a; Abdul-Aziz et al. 1995; Montes 
et al. 1996b, 1997; Montes & Ramsey (1998); Abranin et 
al. 1998; Berdyugina et al. 1998) and without obvious line 



asymmetries in other cases (Hawley & Pettersen 1991). The 
authors ruled out the possibility of pressure (Stark) effects 
as the cause and conclude that the profiles are best explained 
if these broad components and asymmetries are attributed 
to plasma turbulence or mass motions in the flare. In so- 
lar flares, most frequently, a red asymmetry is observed in 
chromospheric lines and interpreted as the result of chro- 
mospheric downward condensations (CDC) (Canfield et al. 
1990 and references therein). However, evidence of a blue 
asymmetry has been also reported (Heinzel et al. 1994) and 
blue and red asymmetries are observed simultaneously at 
different positions in a flaring region (Canfield et al. 1990), 
or at the same position but different times (Ji et al. 1994). 
Recent line profile calculations (Can, Rieger & Fang 1993; 
Heinzel et al. 1994; Ding & Fang 1997) show that a CDC 
can explain both the blue and red asymmetries. On the 
other hand, in stellar flares evidence of mass motions have 
also been reported, in particular a large enhancement in the 
far blue wings of Balmer lines during the impulsive phase 
of a stellar flare was interpreted as a high velocity mass 
ejection (Houdebine, Foing & Rodono 1990), or high ve- 
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locity chromospheric evaporation (Gunn et al. 1994b), and 
red asymetries in the wings of Balmer lines are reported by 
Houdebine et al. (1993) as evidence of CDC. Broadening 
and shifts in UV lines during stellar flares, also interpreted 
as mass motions, have been reported by Simon, Linsky & 
Schiffer (1980), and Linsky et al. (1989). Thus another pos- 
sible explanation of the broad components and asymmetries 
we observed is mass motions in the flare - perhaps an explo- 
sive ejection/eruption (blueshift) followed by flows down the 
flaring loops (redshift). Still, since a CDC can also explain 
both the blue and red asymmetries observed in this stellar 
flare, it remains a distinct possibility as well. 



mass motions in the flare. Similar blueshifts may be seen in 
the C IV line during the impulsive phase of the flare. 

Ultraviolet TR lines are enhanced by a factor of > 20, 
chromospheric lines by a factor of > 4 and the far UV con- 
tinuum by a factor of > 10 x over background in our UV 
spectrum just after optical flare maximum. The continuua 
in our flare-affected UV spectra generally agree with a Si I 
recombination model. 

We estimate the energy released during the flare for all 
the optical chromospheric lines at ~ 5.7 10 33 ergs, indicating 
that this flare in LQ Hya is more energetic in line emission 
than an average flare on a dMe star. 



4 CONCLUSIONS 

We have detected a strong flare on LQ Hya in high reso- 
lution optical spectra (4800 to 7000 A) and IUE SWP ob- 
servations. Such a strong flare is unusual in an early K star 
like LQ Hya. The flare started on 1993 December 22 be- 
tween 2:42 UT (quiescent spectrum) and 4:07 UT (end of 
first enhanced UV spectrum). UV data suggest the impul- 
sive phase began well before the first optical sign at 04:36 
UT. The optical chromospheric lines reached their maxi- 
mum intensity «55 min later, by which time the continuum 
enhancement had sharply decreased. Thereafter, the optical 
line emission slowly decreased in a gradual phase that lasted 
at least until the last observation (07:29 UT). Quiescent C IV 
flux levels were not recovered after ~4 h UT on the following 
day (though a second flare or rotation of the flaring region 
beyond the limb may have affected the results). 

We detected an optical continuum enhancement that 
increased toward the blue (oc A -1,35 ) and reached a maxi- 
mum (36%) during the impulsive phase. The UV continuum 
was enhanced by at least ~ 10 x at 1500 A and 2800 A. 

We analyse the lines by subtracting a quiescent LQ Hya 
spectrum or that of a modified inactive reference star. The 
excess Ha and H/3 emission equivalent widths, in the ob- 
served - reference spectra, increase by a factor of 2.7 and 
5.8 at maximum, respectively, over the quiescent level. The 
Hq/H/3 Balmer decrement is shallower during the flare, 
changing from 3.15 at the quiescent state to 1.46 at the 
impulsive and maximum phases of the flare. We observe the 
He I D3 line, a well known diagnostic of flares, going into 
emission during the flare, reaching an EW = 0.346 A at the 
maximum. We also observe excess emission in He I lines at 
4921.9, 5015.7, and 6678.1 A in the (O-Q) spectra, and in 
other metal lines such as the Na I Di and D2, the Mg I b 
triplet and several Fe I and Fe 11 lines. 

The more intense lines, Ha, H/3, He I D3 and He I A6678, 
exhibit very broad wings in the subtracted profiles. Their 
profile are hence not well matched with a single component 
model; we have used a two Gaussian fit (narrow and broad) 
to characterize them. For all these lines the contribution of 
the broad component to the total EW of the line and their 
FWHM reach a maximum at the impulsive phase and then 
decrease. Moreover, the line profiles are asymmetric, appear- 
ing the broad component blue-shifted in the impulsive phase 
and red-shifted during the gradual phase, with AA between 
the components increasing with time. Mg 11 profiles respond 
similarly. These broad components and asymmetries can be 
attributed to plasma turbulence or to upward and downward 
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